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Abstract. - We present a general sufficiency condition for the presence of multipartite entangle- 
ment in thermal states stemming from the ground state entanglement. The condition is written 
in terms of the ground state entanglement and the partition function and it gives transition tem- 
peratures below which entanglement is guaranteed to survive. It is flexible and can be easily 
adapted to consider entanglement for different splittings, as well as be weakened to allow easier 
calculations by approximations. Examples where the condition is calculated are given. These 
examples allow us to characterize a minimum gapping behavior for the survival of entanglement 
in the thermodynamic limit. Further, the same technique can be used to find noise thresholds in 
the generation of useful resource states for one-way quantum computing. 



In recent years there has been much effort to investi- 
gate the role of entanglement in general physics problems. 
Entanglement is known to be a key resource for quan- 
' turn information, essential for faithful teleportation and 
allowing an absolute secure key distribution among other 
things [1] , and the study of entanglement has been mainly 
from this perspective. However, entanglement is also a 
key foundational issue in quantum mechanics, and has re- 
cently been associated to various phenomena in different 
areas of physics, for example Hawking radiation in cosmol- 
ogy [2] , symmetry breaking in high energy physics [3] and 
in particular, to many areas of condensed matter physics 
such as critical phenomena [4]. In addition, entanglement 
theory has also been helpful in finding the ground state 
for difficult many-body systems [5]. 

All of these results and connections are very intrigu- 
ing, and lead us to ask when and where else entanglement 
exists, and what is its role in the associated phenomena. 
To address these issues, in the first instance, it would be 
very useful to have a general, easy test to see if a sys- 
tem contains entanglement. There are several diflSculties 
in this, especially in many-body physics. First, in nature, 
systems are in thermal states, and calculating the density 



matrix involves the diagonalisation of large Hamiltonians 
which, in general, proves impossible. Second, given this 
state density matrix, it is difficult to test if it is entangled 
or not. 

One situation that helps simplify the problem is given by 
systems where the ground state is highly entangled. This 
is often the case for symmetric or interacting many-body 
systems. There, at low enough temperatures, the prop- 
erties are governed by the ground state, and the system 
is entangled. Along these lines, for example, in [6,7], by 
taking the minimum expectation of the energy allowed for 
separable states, conditions for entanglement are found on 
the average energy and associated thermodynamic quanti- 
ties, so observing energy below this minimum value means 
the system is entangled. 

In this paper we present an explicit connection between 
the ground state entanglement and the entanglement of 
the thermal state. We give a condition for the existence of 
entanglement based on minimum knowledge of the ground 
state and statistical properties of the full state (the parti- 
tion function) , so that we do not need to calculate the full 
density matrix. This is a kind of coarse grain approach 
to the existence of entanglement. Weaker, more easily 
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calculable versions using approximations to these quanti- 
ties follow easily. The methods can further be adapted 
to consider different types of entanglement (from the full 
multipartite entanglement to bipartite entanglement), and 
also give ever more precise entanglement conditions, up to 
the exact case. We present several examples which in- 
clude extreme (good and bad) cases, allowing us to ob- 
serve a minimum requirement on the gapping behaviour 
of the Hamiltonian if we are to see entanglement from the 
ground state survive in the thermodynamic limit. 



Consider a general system in the thermal state 

PT = e ' " /Z, 



where H is the Hamiltonian, T is the temperature, Z is the 
partition function and fcs is the Boltzmann constant. The 
thermal states can be expanded in their diagonal basis, the 
energy eigenbasis, as 

-E,/kBT 

E z 



^-Ea/ksT 

PT = ^ |eo)(eo| 



where |eo) is the lowest energy eigenstate, which is often 
entangled and je^), j > 0, are the (possibly entangled) ex- 
cited states having energies Ej > . For any temperature 
T ^ the thermal state is mixed except for T — when 
Pt becomes identical to the pure ground state (the lowest 
energy eigenstate |eo)). It is apparent that if the lowest 
energy eigenstate, |eo), is entangled, then the system will 
be entangled at T = and it seems reasonable that it will 
retain at least partial entanglement for small, finite tem- 
peratures. The immediate question is, when T grows, and 
a portion of excited states gets mixed in, is the thermal 
state for such finite T still entangled? And what is the 
critical temperature for the survival of that ground state 
entanglement? To answer this question we pick an entan- 
glement measure that suits the problem best: The global 
robustness of entanglement [8] measures the worst degra- 
dation of entanglement through mixing and allows there- 
fore to establish a minimal temperature range in which 
the ground state entanglement survives within the ther- 
mal mixture. 

Definition 1 The global robustness of entanglement, 
R{(j), is defined for a state, a , as the minimum amount, t, 
of arbitrary noise, t, that needs to be mixed to a to make 
the state separable, i.e. 



Ria) 



min t, 

w(T,t) 



(3) 



(4) 



such that there exists a state t satisfying 

Lu{T,t) := ^((T + tr) e S, 

where S is the set of separable states with respect to an 
arbitrary partition of the system into subsystems. If a is 
separable with respect to that partitioning, then no noise 
needs to be added and the robustness is zero, R{(j) = 0. 
Conversely, R^cr) > if and only if a is entangled. 




Fig. 1: The probability of the lowest energy eigenstate, 
po = - — — ! is plotted against temperature T for a typ- 
ical thermal state. The point where the curve crosses the 



(1) value 



defines the transition temperature Tirana, below 



l+fi(l'=o>) 

which entanglement is guaranteed 



Taking a thermal state pT, expanded as in Eq. (2), as 
a candidate for cu in Eq. (4) one can infer how much ther- 
mal add-mixing is at least allowed until all entanglement 
could vanish. This can be achieved by identifying, for in- 
stance, the lowest energy eigenstate |eo) with the state a 
and the higher energy eigenstates as thermal noise. In 
this scenario the pre-factor of a is just the probability of 
the lowest energy eigenstate within the thermal mixture, 
i-S- = ~ — °z ^ • This argument leads immediately 
to a sufficient criterion for the presence of entanglement 
at finite temperature, which is stated in Lemma 1 and 
illustrated in Fig. 1. 

Lemma 1 A thermal state px as defined in Eq. (1), with 
partition function, Z , and whose lowest energy eigenstate, 
|eo), has energy Eq and global robustness of entanglement 
i?(|eo)), must be entangled if 



-Eo/ksT 



> 



1 



l + i?(|eo)) 



(5) 



holds. 



The Lemma follows directly from the definition of the 
global robustness of entanglement and is not surprising 
with respect to the mathematics. However, it is worth 
discussing the physical aspects of the use of this adapted 
measure for thermal states. In particular it is interesting 
to note that relation (5) combines some statistical prop- 
erties of the thermal state, such as the partition function 
and the ground state energy, with the entanglement con- 
tent of the ground state. The fulfilment of condition (5) 
detects or witnesses the existence of entanglement in the 
thermal mixture, even though entanglement can also exist 
when Eq. (5) is not met. Despite this fuzzy transition, it 
is interesting to define a transition temperature Ttrans by 
setting equality in (5), see Fig. 1, which tells when entan- 
glement must necessarily occur in the system. To make a 
first guess whether a thermal state is entangled and give a 
minimal transition temperature only the three ingredients 
stated above are needed - as opposed to the full knowledge 
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of all energy eigenstates or the entanglement properties of 
higher excited energy states. Of course, knowing these 
could lead to more accurate bounds on the temperature 
range where entanglement is present. The condition itself 
however is only tight, when the thermal mixture degrades 
the entanglement of the ground state in the most efficient 
way. Indeed entanglement of the thermal state may go up 
as the temperature increases, however, this would not be 
due to the ground state's entanglement, and it is typically 
not so, except in very special examples. 

Obervations: Let us make a few observations about 
the validity and power of condition (5). Firstly, it can be 
simply adapted to also consider the situations where the 
main entanglement contribution does not come from the 
ground state, but some excited state \ej), or even contribu- 
tions from several states. This is done by replacing i?(|eo)) 
with the robustness of the state in question, R{\ej)), in the 
right-hand-side (RHS) of Eq. (5) and the state's popula- 

tion pj = - — on the left-hand-side (LHS). Similarly 

we can adopt degenerate ground states by replacing the 
state I Co) by the maximally mixed state across the degen- 
erate space. We could also, for example, take the first few 
excited states as the main contributor, right up to consid- 
ering the whole state i.e. including all high energy levels. 
Of course, the price to pay is that the robustness of en- 
tanglement will become more difficult or even impossible 
to calculate. 

Secondly, by considering the entanglement of the ground 
state with respect to different partitions into subsystems 
the condition can witness different types of entanglement. 
A system of many parties can be re-arranged by grouping 
some of the parties together, for instance as a bipartite 
grouping, and the global robustness of entanglement for 
such grouping could be applied in Eq. (5) as well as the 
robustness of the full multi-partite partitioned system. Of 
course the set of totally separable states with respect to all 
parties is a subset of the set of bipartite separable states 
and therefore the robustness decreases the fewer parties 
we distinguish, R{p) > Rsiip)- The LHS of Eq. (5) is a 
monotonically decreasing function of temperature so that 
higher ground state entanglement implies a higher transi- 
tion temperature. Therefore, the corresponding transition 
temperature for a bipartite grouping, T^^^^, must be less 
or equal then the transition temperature for full multi- 
partite entanglement, T^^^^^ < Ttrans- This is similar to 
the situation in [7], where the transition temperature in- 
creases with the number of partitions. Fortunately, Rsiip) 
is known for pure states [8] (for any bipartite cut), so it is 
always possible to calculate a bound to the RHS of Eq. (5) 
by considering only bipartite cuts for a given ground state. 

Thirdly observe, that when the partition function and 
the robustness of entanglement is replaced by upper and 
lower bounds to those quantities, respectively, Eq. (5) is 
still satisfied and presents a valid entanglement condition. 
Thus, we may use approximation methods for all calcula- 
tions, providing they bound in the correct direction. An 



example of such a lower bound to the global robustness of 
entanglement is given in [9]. For any pure state, jV'), 

H-i?(|V')) > 2^«(l'^» > 2^«(l'''», (6) 

where Ef{{\ij:)) is the relative entropy of entanglement [10] 
and Eq{\%1))) is the geometric measure of entanglement 
[11]. These distance-like entanglement measures allow an 
intuitive interpretation of the entanglement condition (5). 
One can understand Eii{\^)) and Eq{\iI))) as the minimum 
"distance" £>([?/') [[w) to the closest separable state lo with 
respect to relative entropy and the geometric overlap re- 
spectively. Thus, ii D{\e^)\\pT) < D{\eo)\\Lo) = Ei\eo)) 
then Pt is entangled (illustrated in Fig. 2). 




Fig. 2: The figure shows the set of all states and the sub-set 
of separable states (shaded area). The distance between the 
ground state (eo) and the closest separable state to gives the 
entanglement of the ground state E{\eo)), illustrated by the 
dotted line. The solid line is the line of thermal states pr for 
increasing temperature from T = (po ~ |eo)(eo|) to poo = ^1 
for T —* OQ. If the distance from |eo) to pr is smaller than 
distance from |eo) to uj, then pT must be entangled. 

For the thermal state, both the relative entropy and 
the geometric measure can be easily calculated giving 
D{\eo)\\pT) = -log(po), where po = e-^°^''"'^/Z is the 
ground state population. This relation reproduces (5), 

Po > 2-°(l"">ll"), (7) 

when i?(|eo)) is replaced with the bounds in (6). The 
statistical interpretation of the relative entropy [12] on 
the RHS of (7) offers another, statistical, interpretation of 
Eq. (5) and Eq. (7): Entanglement is present if the prob- 
ability of the system being in the ground state is greater 
than the (asymptotic) probability that [eo) would be mis- 
taken for its closest separable state uj. 

Example 1: We now consider three illustrative exam- 
ples for which we derive transition temperatures for entan- 
glement. Our first example is a spin dimer, that is, two 
spin-1/2 particles interacting via the Heisenberg Hamil- 
tonian with coupling J, exposed to an external magnetic 
field B, 

H ^ 5(4 + 4) + j(44 + 44 + 44), (s) 

where cr^x/Y/z Pauli matrices for the ith spin and 

B,J > 0. This example allows us to compare our condi- 
tion (5) to the exact case which was calculated in [13]. 
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Diagonalising the Hamiltonian leads to the four eigen- 
values J+B, J, J—B, —3 J from which we can immediately 
calculate the LHS of (5). We have Z = e-("'+'S)/'=«^ + 

e-J/knT + ^-iJ-B)/ksT ^ g3J/fc«T for B < AJ the 

ground state is the singlet with eigenvalue — 3J. The ro- 
bustness of the singlet is known, i?(|eo)) = 1 [8]. So, for 
B < 4 J we have the sufficient condition for entanglement 



-AJ/ksT^^B/kBT 



-B/ksT 



1)<1. 



(9) 



In particular, when _B = 0, the excited states are all degen- 
erate with eigenvalue J, and the entanglement condition 
can be rewritten to give a transition temperature below 
which the system is guaranteed to be entangled 

4J , , 

T < -— - = Ttrans. 10 

Kb In 3 

This remarkably coincides with the exact entanglement re- 
sult in [13]. For non-zero B the excited states lose their 
degeneracy and bound (9) gives a lower transition tem- 
perature than the exact result. This happens because our 
inherent approximation that the thermal mixture kills off 
the entanglement optimally, becomes incorrect. 

When B > 4 J, the singlet ceases to be the ground state, 
and is instead superseded by the separable state |00) with 
energy J — B. The system thus undergoes a quantum 
phase transition. Expression (9) still gives a valid suffi- 
cient condition for entanglement, but now the main con- 
tribution comes from the first excited state (in this case 
the singlet), rather than the ground state (which is no 
longer entangled). However it is no longer satisfied for 
any temperature, and so does not see the presence of en- 
tanglement, even though we know it exists from [13]. 

Example 2: The second example is a simplistic very 
general model for a class of many-body Hamiltonians, 
which allows us to examine extreme cases, and make in- 
teresting general statements. We only specify the energy 
spectrum (giving the LHS of (5)) and the ground state, 
common in many-body physics (whose robustness gives us 
the RHS of (5)). This toy model illustrates the flexibility 
we have in considering these quantities separately. 

The ground states we consider arc of the form 



\Sin,k)) 



1 



PERM 



(11) 



where the sum is over all permutations of position. These 
states describe n systems with k excitations symmetrised 
over all permutations. We could, for example, imagine n 
sites, which can either be empty or occupied by k iden- 
tical (fermion) particles. These are typical ground states 
in solid state physics, for example in some Hubbard and 
related models. They also arise as Dicke states when n 
atoms interact with a single electromagnetic mode [14]. 

It can be shown [15] that the symmetry of these states 
means that the global robustness of entanglement is re- 
lated to the relative entropy of entanglement {Er) as 

1 + R{\Sin, k))) = 2^^«(l^^("^'=)». (12) 



The relative entropy of entanglement for the states 
\S{n, k)) is known [14, 16] and gives 



1 



R{\S{n,k)))^l^Q 



n — k 



(13) 



For simplicity we now only consider the case where k = 
n/2. This choice of k maximises the entanglement and for 
large n we have 



1 -t- Ri\Sin, n/2))) = 2^«(l^("''=)» = y/Ti. 



(14) 



For the LHS of (5) we consider Hamiltonians of the form 



D-l 



F-So|eo)(eo|+ ^(^o + m"A)|e™)(e, 



(15) 



where D is the number of energy levels (the dimension 
of the space), with spacing parameter < a < 1 and A 
is the energy gap between the ground state and the first 
excited state. This gives a range of energy spectra from 
the case where all excited states are degenerate {a ~ 0) to 
the case where they are equally spaced (a = 1). Although 
this model is rather simplified and hence somewhat arti- 
ficial, as mentioned, its simplicity allows us to examine 
extreme cases, through which we arc able to make general 
statements of interest true for all systems. 

For simplicity here, we will use the relative entropy of 
entanglement to consider transition temperatures. As in 
the above entanglement calculation, it will not effect any 
of our conclusions to use this weaker condition. 

The thermodynamical partition function is 



D-l 



,-m°A/fcBT 



(16) 



For a = 0, the degenerate case, we can easily evaluate the 
sum in (16), giving the condition for entanglement 



T < 



A 



1 



iB In 



D-l 
2Eh(I = o>)_i 



(17) 



In fact Zq is an upper bound to any partition function. 
Hence Tq in (17) is a lower bound to all possible threshold 
temperatures, written only in terms of the gap to the first 
excited state A and the entanglement of the ground state. 

However, in the thermodynamic limit, this lower bound 
will almost certainly tend to zero. We can see this by not- 
ing that the entanglement of any state is trivially upper 
bounded by logD. This is found by considering the dis- 
tance to the maximally mixed state. In general this value 
cannot be achieved [17], and so E < logD, and at best 
we can hope for the dimension over entanglement term to 
give a constant in scaling. More usual (in fact as far as 
the authors know, there are no counter examples) , will be 
as in the example state (14) where the entanglement is so 
low that this term will cause the temperature to tend to 
zero in the limit. 
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Of course we must also consider the scaling of A, but 
again, at best this is likely to be constant, and at worst 
(and more usually) it will also limit to zero in the thermo- 
dynamic limit. We can conclude that in this a ~ "worst 
case" spectrum, even with a highly entangled ground 
state, the entanglement cannot survive any non-zero tem- 
perature in the thcrmodynamical limit. 



For a = 1 we get Zi 



-Eo/ksT 



D-l 



the limit of D — > oo we find the entanglement condition 



T < 



A 



1 



iB In 



2^h(I'=°)) 

2Efj(l<!o>)_l 



Tl. 



(18) 



When we take a low temperature approximation of Zi the 
transition temperature Ti scales roughly with 2^"^!^°^^ 



(19) 



So, in the limit of large n, for fixed A, we have entan- 
glement for arbitrarily high temperatures. This spectrum 
is in some sense a best case extreme, and we would not 
expect such a nice situation in real systems. However it 
does allow us to see that if we wish to see entanglement 
in the thcrmodynamical limit, the scaling of A should not 
be too bad compared to the scaling of the entanglement. 

Thus, to see entanglement, even for the best case energy 
spectrum, we require an energy gap between the ground 
and first exited state that scales at least of the order 



A > 



1 



(20) 



Any less would mean that the transition temperature 
tends to zero in the thermodynamic limit. 

If we take the example state entanglement (14), we see 
the transition temperature 



Kb 



(21) 



Now, the way A scales with n is very dependant on the 
system. To see entanglement we would want a system with 
scaling A > l/y^. 

We now consider the intermediate case for completeness. 
For a 7^ and in the limit D — > oo we can approximate 
the sum in (16) by an integral which has the same form 
as Riemann's gamma function F giving 



-Ea/ksT 



r(^) (kBT 



l/a 



(22) 



Thus, with (14) we get the condition for entanglement 

= T,. (23) 



^ kB 



F i 



For a = \ this reduces to the low temperature approxima- 
tion (21). Similarly, for all a > the transition temper- 
ature increases with n. Indeed for any ground state with 



entanglement scaling with n a similar statement holds. 
Again we see, that for our technique to see the survival of 
entanglement from the ground state, we would need some 
limits on the scaling of the gap A. This is in analogy 
to many situations in condensed matter physics, where to 
witness interesting phenomena we need a gapped system. 

As mentioned these examples are somewhat artificial, 
and realistic Hamiltonians will lie somewhere in between 
the a = and a ~ 1 cases. However, with the fact that 
the gap plays such an important role in both extremes 
(and for in between a), and the fact that as soon as we 
step out of the a = case we have non-zero threshold 
temperature, even when the gap scales to zero in the limit 
(albeit conditioned that it does not do so too fast), we 
may conjecture that for any realistic Hamilitonian, with 
fixed gap A, any scaling of the entanglement of the ground 
state with D implies a non-zero threshold temperature. 

Example 3: Our third and final example is that of 
the so called stabiliser Hamiltonian, whose ground state 
is an important resource for quantum information, most 
notably for one-way quantum computing [18]. 



H 



K, 



(24) 



where the Ki are the so called generators of the stabilizer 
group, and N(i) are the neighbours of i in a given lattice 
(or graph). Although this Hamiltonian is not a naturally 
occuring one in many-body physics, there are proposals 
for its implementation. Further, its thermal state is also 
identical to the mixed state induced by several realistic 
noise models for the generation of cluster states in various 
schemes, for example using optical lattices, (see e.g. [19]). 

The ground state of this Hamiltonian is the graph state 
of the associated lattice (or graph) . The excited states are 
achieved by making local al, flips, acting as an excitation 
operator. For all lattices, the ground state has energy 
~nB, and the ith excited state has energy B{—n+2i) with 
degeneracy (") . The partition function is easily found as 



Z 



iB/kT 



1 



„2B/fcT 



g2B/fcT 



(25) 



Our condition for the survival of entanglement then reads 
(again taking relative entropy of entanglement) 



T < 



2B 



(26) 



fcln(2^«(l'=o))/" - 1) 

This can also be translated into conditions on allowed 
noise when directly constructing the graph state. For ex- 
ample, suppose a local cr* error occurs with equal proba- 
bility P on any site on such a graph state, this noise state 
can easily be seen to corresponding to the thermal state of 
(24), with P = i_^_^2B/kT ■ Hence this translates to a con- 
dition on tolerable flip error threshold, below which the 
state is still entangled, of 



P < 1 - 2~-^«(l^">)/" = Ptr 



(27) 
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For many graph states, including the ID, 2D, and 3D 
cluster states, the entanglement is calculated to be Er = 
n/2 [20], giving threshold temperature and noise level of 



2B 



/cln(V2- 1) 



Ptrans — 1 1 



(28) 



Incredibly this coincides exactly with the temperature 
below which entanglement can be distilled [19]. Even 
more, in the same paper Kay et al. also show that above 
this temperature entanglement cannot be distilled. This is 
very encouraging, since by completely independant meth- 
ods we have arrived at the same transition values. Thus, 
here our condition not only implies the existence of entan- 
glement it also turns out that the detected entanglement 
is of the useful, distillable, kind. Whether this is a gen- 
eral property of the criterion is presently unclear nonethe- 
less an intriguing possibility. Finally, our approach will 
work for other lattice shapes and perhaps allow extensions 
where the methods of [19] fail. This is beyond the scope 
of this work, and we rather leave it for future study. 

In conclusion, in this paper we have presented a gen- 
eral condition for entanglement in thermal states, based 
on the entanglement of the ground state. The condition 
separates the statistical physics part from the entangle- 
ment part, so that results may be used from either field 
together to give transition temperatures below which en- 
tanglement is guaranteed. The flexibility of the condition 
allows for approximations, allowing easier calculation, and 
for different kinds of entanglement to be studied. Further, 
by considering extreme examples we have seen that even 
in the best case (in terms of the energy spectrum, thus Z) , 
we need some gapping if entanglement is to survive in the 
thermodynamic limit, as characterized by equation (20). 
With the final example we see that the approach can be 
used, not only in many-body physics, but also for simu- 
lating noise models for prcperation of resource states for 
quantum information, to find tolerable noise thresholds. 

Through Z and the equations of state, the main condi- 
tion (5) can be rewritten in terms of other thermodynamic 
quantities as in [7] , allowing the prospect of a variety of ex- 
perimentally accessible quantities detecting entanglement. 
These results should provide a very useful tool for investi- 
gating entanglement in many statistical physics systems. 

Several natural questions arise from these investiga- 
tions. Firstly, with respect to quantum information, we 
can ask, given that there are systems which can contain 
entanglement in the thermodynamic limit, how might we 
access this resource? This is the topic of several recent ef- 
forts, for example by scattering [21], and is the topic of on- 
going investigations. Secondly, with respect to condensed 
matter issues, we can ask how these relate to critical phe- 
nomena. In this respect it is intriguing that our results 
indicate the requirement for certain gapping to observe 
the survival of entanglement. Though this and the gen- 
eral connection to entanglement remain large open ques- 
tions, as mentioned, with the connection to other thermo- 



dynamic quantities, we hope that these connections can 
be pushed further and our understanding increased. 
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